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Abstract 

The nonlinear evolution of a pair of initially linear oblique waves in a 
high-Reynolds-number Stokes layer is studied. Attention is focused on times 
when disturbances of amplitude e have O^^R) growth rates, where R is 
the Reynolds number. The development of a pair of oblique waves is then 
controlled by nonlinear critical-layer effects (Goldstein & Choi, 1989). Vis- 
cous effects are included by studying the distinguished scaling e = 0(R -1 ). 
This leads to a complicated modification of the kernel function in the integro- 
differential amplitude equation. When viscosity is not too large, solutions to 
the amplitude equation develop a finite-time singularity, indicating that an 
explosive growth can be induced by nonlinear effects; we suggest that such 
explosive growth can lead to the bursts observed in experiments. Increasing 
the importance of viscosity generally delays the occurrence of the finite-time 
singularity, and sufficiently large viscosity may lead to the disturbance de- 
caying exponentially. For the special case when the streamwise and spanwise 
wavenumbers are equal, the solution can evolve into a periodic oscillation. A 
link between the unsteady critical-layer approach to high-Reynolds-number 
flow instability, and the wave-vortex approach of Hall & Smith (1991), is 
identified. 


1 Introduction 


The flow generated above a sinusoidally oscillating plane wall is one of the simplest 
exact unsteady solutions of the Navier-Stokes equations. The instability of such 
Stokes layers has been studied as a paradigm of the instability of unsteady (periodic) 
flows. Although the flow is unidirectional, a conventional normal-mode approach 
to the linear stability of the flow is not possible due to the unsteadiness of the basic 
state. Instead von Kerczek Sc Davis (1974) and Hall (1978) used Floquet theory to 
seek linear disturbances which grow over a complete period. However, they found 
that over a full period the flow was stable at all Reynolds number investigated. At 
the highest end of the range studied, this included Reynolds numbers for which 
instabilities have been observed experimentally (e.g. Merkli &: Thomann, 1975). 

This paradox has been partly resolved by Tromans (1977) and Cowley (1987) 
who argued that at high Reynolds numbers the rapid growth of small high-frequency 
disturbances over part of a period can lead to nonlinear effects preventing the linear 
decay over a whole period (see also Hall, 1983). This idea was developed by Wu 
Sc Cowley (1992) (see also Wu, 1991) for two-dimensional disturbances using the 
unsteady, or non-equilibrium, critical-layer approach of Hickernell (1984), Churilov 
Sc Shukhman (1988), Goldstein Sc Leib (1989) and others. They found that following 
a linear development stage, nonlinear interactions inside the critical layers could 
affect the evolution of disturbances sufficiently to cause the amplitude to ‘blow-up’ 
in a finite time. In this paper we extend our analysis of instability and transition 
in Stokes layers to three-dimensional disturbances consisting of a pair of oblique 
waves. 

The importance of three-dimensionality has long been realised in transition. For 
instance, in the case of the Tollmien-Schlichting instability of boundary layers, dis- 
turbances are predominantly two-dimensional only in carefully controlled exper- 
iments, e.g. Schubauer Sc Skramstad (1947), Nishioka et al. (1975). Even then 
two-dimensional disturbances dominate only in the early stages of transition, with 
three-dimensional disturbances growing to significance downstream, e.g. Klebanoff 
et al. (1962), Kachanov Sc Levchenko (1984), Kachanov (1987), Saric Sc Thomas 
(1984). For the Stokes layer, three-dimensionality appears to be even more impor- 
tant because of the apparent difficulty of introducing, and then observing, artificial 
two-dimensional disturbances; indeed to the best of our knowledge there are as yet 
no experiments on Stokes layers with controlled two-dimensional disturbances. 

In order to analyse the three-dimensional instability of Stokes layers, there are a 
number of theoretical approaches available. Since no finite critical Reynolds number 
has yet been identified from a linear Floquet theory analysis, a weakly nonlinear 
expansion about the critical Reynolds number in the spirit of Seminara Sc Hall 
(1977) is not possible. Instead we assume that the Reynolds number, R , of the flow 
is large so that the frequency of the Stokes layer, u>, is much smaller than a typical 
0(ujR) frequency of the instability waves. Under these conditions, linear instability 
waves are quasi-steady and satisfy Rayleigh’s equation (Tromans 1977). Figure 1 is 
a graph of the neutral curves of a Stokes layer plotted as parametric functions of 
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time. This is for a flow where the velocity of the boundary y* = 0 is ( Uo cosu/t*, 0, 0), 
(x*,y*,z*) and t* are dimensional Cartesian coordinates and time respectively, 



are the Reynolds number and Stokes layer thickness respectively, and u is the vis- 
cosity. The streamwise and spanwise wavenubers that are nondimensionalised by 
S*~ 2 are denoted by a and ±/3. At any time an ‘infinite’ number of Rayleigh modes 
exist, but the most rapidly growing modes can for the most part be found at those 
times and for those wavenumbers that lie beneath the solid curve A in figure 1; 
see Cowley (1987) for further details including the significance of the mode crossing 
point marked by X. 

In this paper we will not be concerned with how Rayleigh modes are excited in 
the Stokes layer, i.e. the receptivity problem. Instead we follow other authors and 
assume that as a result of background disturbances Rayleigh modes are introduced 
into the flow. For instance, suppose that a pair of oblique neutral modes are intro- 
duced into the Stokes layer at a time and with a wavenumber given by the left-hand 
branch of one of the curves A or B. As the Stokes-layer’s mean profile slowly evolves 
this mode will begin to grow exponentially fast. A number of possibilities then arise: 


1. First suppose that although small, the initial amplitude is such that nonlinear 
effects come in almost immediately near the left-hand branch of either curve 
A or B. For definiteness we assume that the disturbance is introduced at the 
neutral time r = r 0 with an amplitude e<j. At a time r, where At — t — t 0 <C 1, 
a simple linear quasi-steady theory gives the growth rate as & 0 RAt and the 
wave amplitude, e, as 

c ~ €o exp ^-d- 0 RAr 2 j , 

where &o is a known constant (Cowley, 1987). Inviscid arguments similar to 
those of Goldstein & Choi (1989) then show that nonlinear effects become 
important in an unsteady critical layer when the growth rate is 0(e) R). Of 
course this scaling only applies if the flow becomes nonlinear while a critical 
layer exists, i.e. if At < 1, or equivalently if 

2 log Cq 1 <C croR ■ 


In addition, the initial amplitude e 0 must not be too large, otherwise viscous 
effects must be included in the critical layer. In particular, if e 0 is tuned so 
that nonlinear effects come into play when e = 0(R~ 1 ) and At = 0(R ~ ?), 
i.e. if 

\o Z e^ = 0(R^), 


then the critical layer is unsteady and viscous. This is the scaling that will 
be studied here. The recent excellent comparison of Hultgren (1992) between 
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experiments and a theory based on similar assumptions suggests that theoret- 
ical results so obtained may have a wider range of validity than might first be 
thought in view of our specific choice of distinguished scaling. We delay until 
later a discussion of the Very viscous* case that occurs when the flow becomes 
nonlinear for A r <C i.e. when log £q 1 < Ri (see §4.2 & §4.4). 

2. Next we suppose that the initial amplitude of the disturbance is sufficiently 
small so that the modes evolve linearly until they approach the right-hand 
branch of curve A, 1 say at r = r 0 again. Then, as above, nonlinear ef- 
fects can be accounted for by an unsteady, viscous, critical-layer analysis if 
fH < | Ar| <C 1 (see §4.4 for a brief discussion of the very viscous case 
Ar = 0(12-2)). 

3. The third alternative is that for disturbances of a given wavenumber, nonlinear 
effects become significant when the growth rate of the modes is 0(R) y i.e. at 
a time far away from the neutral curves A and B. A fully nonlinear theory 
then seems necessary, and we do not consider this possibility further here. 
However, we note that there are some asymptotic theories that modify the 
mean flow, so effectively changing the position of the neutral curves, e.g. the 
wave-vortex interaction theory of Hall & Smith (1991). The latter theory 
assumes that there is a weakly nonlinear, high-frequency, neutral wave (e.g. 
a Rayleigh neutral mode) whose evolution on a ‘slow’ timescale nonlinearly 
modifies the basic state by an order one amount. However, in order for this 
interesting theory to be applicable, it is necessary for the flow to be able to 
evolve to the weakly nonlinear neutral wave in the first place. It turns out that 
the viscous generalisation of Goldstein & Choi’s (1989, 1990) work is related to 
this question, and provides a link between two of the current high-Reynolds- 
number stability theories (see §4.4). 

A main concern of this paper will thus be with the nonlinear effects associated 
with unsteady, viscous critical layers. Specifically we shall seek a possible mechanism 
for the bursting phenomenon observed in experiments (e.g. see Merkli & Thomann 
(1975), Hino et ah (1976, 1983)). We shall assume that the disturbance consists of 
a pair oblique waves because the ‘quadratic interaction’ of such waves produces a 
vortex flow, and we note that a significant vortex structure has been observed by 
Hino et ah (1983). 

There have been two main strands in the critical-layer theory of nonlinear flow 
stability (see the reviews by Stewartson (1981) and Maslowe (1986) for overviews 
of critical-layer theory and its applications). One strand has consisted of directly 
seeking nonlinear neutral-wave solutions, but without giving detailed consideration 
to whether the flow could evolve to these equilibrium states, e.g. Benney & Bergeron 

1 The left-hand branch of such modes is on curve A for wavenumbers above point X of figure 1, 
and on curve B otherwise (Cowley, 1987). 
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(1969), Smith &; Bodonyi (1982a, fc), Bodonyi et ah (1983), Gajjar h Cole (1989) 
and Gajjar (1990) (see Gajjar & Smith (1985) for an analysis of the stability of 
such neutral states). In the second strand, the evolution of a linear wave is fol- 
lowed into the nonlinear regime by introducing a slow time (temporal instability) 
or equivalently a ‘slow’ lengthscale (spatial instability), e.g. Churilov & Shukhman 
(1987a, 1988), Goldstein, Durbin h Leib (1987). We take the second approach here, 
and as is appropriate for Stokes layers, we adopt a temporal instability viewpoint. 
However, there are close analogies between the temporal instability of Stokes layers 
and the spatial instability of shear layers and (compressible) boundary layers. This 
allows us both to draw on previous work and to extend it (see Appendix B). 

The instability modes we consider have a wavelength comparable with the thick- 
ness of the Stokes layer. In the case of two-dimensional modes, the scalings and 
nature of the nonlinear analysis depend crucially on whether the neutral modes 
have a regular or (logarithmically) singular critical layer. In the former case, which 
is so on most of the left-hand branch of curve A, a strongly nonlinear critical-layer 
analysis is necessary - cf. Churilov & Shukhman (1987a), Goldstein, Durbin & Leib 
(1987), Goldstein & Leib (1988), Leib & Goldstein (1989), Goldstein h Hultgren 
(1989) and Hultgren (1992) (for other examples of this type of critical layer see the 
work of Goldstein & Wundrow (1990) and Shukhman (1989)). However, on the rest 
of curve A, and all of curve B, there are two or more singular critical layers. Close 
to these curves it is found that for a given growth rate, nonlinear effects are felt 
at smaller disturbance amplitudes than if the critical layers had been regular; as 
a result a weakly nonlinear analysis is possible. However, the amplitude equation 
is not a Stuart-Watson-Landau equation unless viscous effects are large; instead an 
integro-differential equation is recovered (Hickernell (1984), Churilov & Shukhman 
(1988), Shukhman (1991), Wu (1991), Wu & Cowley (1992)). 

In the case of three-dimensional disturbances, Benney (1961) observed that an 
nearly neutral Rayleigh mode has a pole singularity in the streamwise velocity at the 
critical layer. At first sight this suggests that a different scaling is necessary in order 
to follow the nonlinear evolution of three-dimensional modes. Nevertheless, in the 
case of a single oblique mode in a (compressible) shear layer, a Squire transformation 
enables the weakly nonlinear problem to be reduced to one in which the scalings 
are essentially those for a two-dimensional singular critical layer (Goldstein &; Leib 
1989, Leib 1991). This is not possible in the case of a pair of oblique modes. 

The nonlinear spatial evolution of a pair of oblique waves in a free shear layer 
has been studied by Goldstein & Choi (1989, 1990) (henceforth referred to as GC). 
They showed that because of the pole singularity, nonlinear effects must be included 
when the linear growth rate is 0(e*R). Our temporal analysis of the development 
of a pair of oblique waves closely follows the spatial analysis of GC, although we 
additionally allow the critical layer(s) to occur away from inflection points. Even 
though such a generalisation leads to a completely different critical-layer structure 
for two-dimensional disturbances, it does not affect the critical layer dynamics in the 
three-dimensional case if the spanwise wavenumber (3 6* (Wu, 1992a). In addition 

we incorporate viscous effects so that the critical layers involved are unsteady and 
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viscous in nature. 

We are thus interested in the evolution of a pair of high-frequency oblique modes 
when nonlinear effects become important near a neutral curve, i.e. either soon after 
the modes become unstable or just before they stabilise. As explained in detail by 
GC, Wu (1991) and Wu &; Cowley (1992), it is appropriate to concentrate on times 
close to 

r = t 0 + e*Ti , 

for some suitable t x = 0(1), i.e. times at which the linear growth rate is 0(e*R). 
Therefore we introduce the time scales 

h = \^ Rt > ( 1 - 1 ) 

and 

t = Rr (1.2) 

to account for the ‘slow’ nonlinear growth/decay of the disturbance, and the ‘fast’ 
carrier wave frequency of the disturbance, respectively. 

The basic flow U evolves on the very slow time scale r, and it turns out to be 
sufficient to express its profile at time r as a Taylor series about the neutral time tq\ 

U(y,r) = U(y,T 0 ) + e>U r (y,T 0 )Ti + ... . 

Hereafter all quantities associated with the basic flow will be evaluated at r 0 unless 
otherwise stated. 

In order to maintain maximum generality, we wish to force viscous diffusion terms 
to appear at leading-order in the critical-layer equations. An elementary balance 
of the unsteady, u T , and viscous, Uyy , terms in the critical layer of width e* (see 
(2.30)), shows that we require 

iT 1 = Ae , (1.3) 

where the parameter A is introduced to reflect the importance of viscosity. It will 
be assumed to be of order one in §2 and §3. The highly viscous case corresponding 
to A being asymptotically large will be discussed in §4.2 and §4.4. 

The overall evolution of a three-dimensioned disturbance is summarised in figure 
2 for the case when the flow goes nonlinear near the right-hand branch of curve 
A. As illustrated, the disturbance is initially linear and grows exponentially until 
its growth rate decreases to 0(e5i?) when t x — 0(1). At this stage, nonlinear 
interactions inside the critical layers control the evolution. We wish to emphasise 
that there are four times scales illustrated in this figure: 

1. The very slow time scale, r, over which the Stokes layer evolves; 

2. The slow time scale, t x> over which the growth rate evolves; 

3. The faster time scale, f 1} over which the disturbance grows; 

4. The fast time scale, i, over which the disturbance oscillates. 
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We note that although our analysis is based on being close to either the left or 
right hand branches of the neutral curve, it is straightforward to generalise it to 
wavenumbers close to the apex of curve A in figure 1 (cf. Hickernell, 1984). 

The paper is organised as follows. In §2 we construct asymptotic perturbation 
expansions in the ‘outer* region away from the critical layers. The limiting forms of 
these solutions near the critical layers are then determined; as usual these contain 
unknown ‘jumps* across the critical layers. A solvability condition is also deduced for 
an inhomogeneous Rayleigh equation. In §3, we analyse the unsteady, viscous and 
weakly nonlinear flow within the critical layers. By matching the inner and outer 
solutions the unknown jumps are evaluated. Then by combining the solvability 
condition with these jumps, we derive the amplitude equation which is a main 
result of this paper. The amplitude equation is studied in §4, both analytically and 
numerically. In particular, a finite-time singularity structure is identified as in GC, 
and confirmed by numerical solution. In addition, exponentially decaying solutions 
are found under certain conditions. The viscous limit is discussed and a link is 
established with the wave-vortex interaction work of Hall & Smith (1991). Finally, 
in §5, we summarise our main results, and discuss the implications of this study. In 
appendix B, we deduce the amplitude equation for free shear layers by combining 
the present results with those of GC. 


2 Outer Expansion 

We take the flow to be described by Cartesian coordinates (z*, 2 /%z*) = £*(x, j/, z), 
where x* is parallel to the direction of oscillation of the plate, y * is normal to the 
plate and z* is the spanwise direction. We non-dimensionalise time with u? _1 , i.e. 
r = ust*, and write the velocity as U 0 (lf, V , W ). Then the basic Stokes-layer solution 
for flow over an oscillating plate is 

(17, V, W) = ( U , V, W) = (cos(r - y)e~\ 0, 0) . 

We denote the perturbed flow by 

(U^W) = {U + u,v y w) . 


2.1 Asymptotic Solutions Near Critical Levels 

Outside the critical layers, the unsteady flow is basically linear and inviscid. It 
is governed, to the order of approximation required in this study, by the inviscid 
equations 

du dv dw _ 

ai + ^ + a7 = 0 ' 


2 R~ 


! du 

dr 


- du dU 

+ u di +v ~si 


dp 
~dx ’ 
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( 2 . 2 ) 


(2.3) 



and 

§ = 5 «fo) + iAMn I (r„) + ... . (2.13) 

For simplicity we have assumed that the two oblique waves are of equal amplitude. In 
principle it is straightforward to extend the analysis to unequal amplitudes; however 
the asymmetry in the amplitudes complicates the algebra, especially in the viscous 
case. Note that as in Wu &: Cowley (1992), the dependence on the slow time scales r 
and T\ is parametric and will not be written out explicitly. The function Vi satisfies 
Rayleigh’s equation 

(U - c)(D 2 - a 2 )n - UyyVr = 0 , (2.14) 

where 

a = (a 2 + ft 2 )* . 
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The boundary conditions are that v x = 0 on y = 0, and v x — ► 0 as y — * oo. 

We let 7 ] = y — where is the j-th critical level at which U = c. Then as 
rj — > ±0, V\ has the following asymptotic behaviour 


v, ~ + tflfo + Pi'#, log |i)|] , 


(2.15) 


where 


1 


4>a = V + ^Virf + 


, and = 1 + qjT} 2 + ... . 

The function v 2 takes the following form 

v 2 = v 2 (y, ti)E cos f3z + v 2 0 ' 2 ^ cos 2 0z + c.c + . . . , 


(2.16) 


where a relatively large longitudinal-vortex component, U 2 ° ,2 ^cos 2 /?z, has had to be 
included in order to match to the inner solution (see §3). The function v 2 is the 
deviation of the eigenfunction from its neutral state, and satisfies the inhomogeneous 
Rayleigh equation 


U 


yy 


_ c)i 

The asymptotic behaviour of v 2 as y — * yi is 


Z OiU yy T Tl A 


U-c 


}vi • (2.17) 


v 2 ~ -bfrj log \t]\ + ( afrj + bfsj) rjlog \r)\ + 
+cf(f> a + df [4> b + pj(j) a log |?7 1] , 


where 


Vi 




Ur 


yy 


Ur, 


11 Ur, 

+ X- 


U 2 

w yy 


2 ' 2 o , u; ’ 




(2.18) 

(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


U 2 

V 

dA - UyUyyy ~ U 

+ (___ !Qt / rn A ) w }. 


( 2 . 22 ) 


Recall that all the basic-flow quantities are evaluated at time To and at the critical 
level y*. The jumps ( a .+ — a"), etc., will be determined by analyzing the critical 
layers in §3. 

From the continuity equation, we can write 


w x = Aw x E sm j3z + c.c . , 
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where w\ satisfies the equation 


dwi 

dy 


+ 


Uy 

U-c 


Wi = —(3vi . 


This has the solution 

u>i = a -1 sin 0 ~ ^i - ^i,yj ■ (2.23) 

where 

sin 0 = (3 /a . 

The velocity U\ has the form 

ui = AuiE cos (3z + u[°' 2 \y, ti) cos 2 f3z + c.c. , (2.24) 

where xii is obtained from the continuity equation as 

y -1 = -(ioty 1 Vi - v 1>y ] sin 2 6 + v liV } . (2.25) 

Note that in order for ui to be able to match with the inner solution (see §3), a 
spanwise-dependent mean flow, u[°’ 2 \y y fi) cos 2/?z, has to be included at leading 
order (see also GC). As will be shown later, this mean flow is driven by a stream- 
wise slip velocity across the critical layer, which itself is generated by a nonlinear 
interaction inside the critical layer. Although the mean flow is large in the sense 
that it has the same magnitude as the fundamental waves, it has no back effect 
on the critical-layer dynamics. Its dependence on t\ can be viewed as forcing the 
longitudinal vortex represented by v[ 0,2 ^ and w^' 2 ^ - 

Similarly, we write the leading-order pressure perturbation as 

Pi = A(ti)piE cos f3z + c.c. , 

where 

pi = ioc~ l cos 6[UyV i — (U — c)v 1(V ] . (2.26) 

As y — > y^, the asymptotic solutions of pi,Ui,Wi become 

pi ~ ia~ 1 U y cosdbf + ... , (2.27) 

Ui ~ — (ia;) -1 sin 2 9bfri~ l + . . . , (2.28) 

wi ~ a -1 sin 6b±rj~ l + . . . . (2.29) 

Note that the singularity in U\ is a simple pole rather than the logarithmic branch 
point characteristic of a two-dimensional (singular) disturbance. As GC observed, it 
is this difference that results in the faster nonlinear evolutionary time scale compared 
with the corresponding two-dimensional case. 
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We now introduce an inner variable: 

Y=\. (2.30) 

e* 

The outer expansions written in terms of this inner variable are then: 

v ~ ebfAE cos (3 z + e* log e^(-bfrj + bfpjAY)E cos (3 z 

+ c *[(-&i r jl°g 1*1 + df ) + M a f Y + bfpjYlo g \Y\)]E cos (3 z 
+ ci log e* [( af rj + bfsj + df P j)Y + ±A Pj bf Y 2 ]E cos f3z 

+ cj [cfY + ( afrj + bfsj + dfpj)Y log \Y\E cos (3z + c.c + . . . , (2.31) 
u ^ e?(— ia) _1 sin 2 6 Ab±Y~ l E cos (3 z + c.c. + • • . , (2.32) 

w ~ e*a -1 sin ObfAY^E cos (3z + c.c + . . . , (2.33) 

p ~ e ia~ l U y cos 0.46* £ cos /3z + c.c +.. . . (2.34) 

2.2 Solvability Condition 

By multiplying both sides of (2.17) by Vi, integrating from 0 to +oo with respect 
to y, and using the asymptotic solutions (2.15) and (2.18), we obtain a solvability 
condition for (2.17): 

^ + Jit,a = - Y. {(*/ c * - b j c j ) - a t - b J a > ) 

-MbfJf - bjij) - (atdf - ajdj )} . (2.35) 

where the sum is over all critical layers, and J\ and J 2 are constants defined by the 
following integrals respectively 


Jl = 

r+00 TJ 

1 Uyy 

Jo (U-c) 2 ldy ’ 

(2.36) 

J 2 = 

/ +00 [ -1- F yyr }v\dy . 

Jo 1 (U~c)V ( U-c) 1 1 y 

(2.37) 


Notably these integrals are singular; the analysis of the critical layers in §3 shows 
that they should be interpreted in the sense of Hadamard. 

After the jumps (at — a^ - ), etc., are determined in the next section, the amplitude 
equation can be derived from (2.35). The nonlinearity is introduced into the am- 
plitude equation through the jumps; thus for the purpose of deriving the amplitude 
equation, we only need to consider those parts of the inner solutions contributing 
to the jumps. This consideration simplifies the algebra to a certain extent. 
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3 Inner Expansion 

Equations (2.31)-(2.34) suggest that the inner expansions within the j-th critical 
layer take the following form 


u 

= e 3 U\ + e 3 U2 + c* U3 + . . . , 

(3.1) 

V 

= + e*V 2 + e$V 3 + ... , 

(3.2) 

w 

= elW l + elW i + eiw, + ... , 

(3.3) 

p 

= P\ + c 3 P2 + e 3 P3 + • • • , 

(3.4) 


where 0(e n loge*) terms have not been explicitly included. This is because as far 
as deriving the amplitude equation is concerned, they are passive in the sense that 
they match onto the outer solutions automatically whenever the solutions at 0(e n ) 
match, i.e. matching at 0(e n loge*) does not yield any additional jump conditions. 
However, these terms must be included if a quantitative comparison between theory 
and experiment is to be made. 

The function V\ satisfies the equation 

^ = 0 , (3.5) 

where 

= S + + < 3 ' 6 > 

The solution which matches the outer expansion is 


V\ = A(ti)E cos (3z + c.c . , (3-7) 

where A = bjA y and 6+ — bj = 6^, i.e. the jump (£>+ — bj) is zero. 

The expansion of the y-momentum equation gives 


sp '-o 

dY ’ 

(3.8) 

and so the appropriate solution is 


Pi = ia~ l U v cos 6AE cos /3z + c.c. . 

(3.9) 

The function Wj satisfies 

T W 9Pl 

LqW\ — 

az 

(3.10) 

We let W\ — W\E sin f3z + c.c., then W satisfies 


L^Wi = iUy sin 6 cos 9 A , 

(3.11) 

where 

Lo ] = dti + nia (U» Y + Urn) X dy2 . 

(3.12) 
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Equation (3.11) can be solved using Fourier transforms to yield the solution 

r( 0 ) 


where 

and 


W, 


W\ — iUy sin 9 cos 9Wq , 

/* "foo 

o (n) = / CAih-Oe-'* ~ ini dt , 

Jo 


n = a(U y Y + U t t x ) , Sl = ^A a 2 U 2 


Similarly, the leading-order streamwise velocity U\ can be written as 

Ui = il\E cos j3z + c.c. . 

It follows from the continuity equation that 


(3.13) 

(3.14) 

(3.15) 



U \ = —Uy sin 2 9Wo°^ . 

(3.16) 

At 0(e*), V 2 satisfies 



TV L V 1 9 f dSu | dS "1 

LM.YY - i,V, + &Y { g x + gz ] , 

(3.17) 

where 




1 d d ^2 ^ 

+ UyrTxY + ~ Urr T l ) + X Qy 2 + ^yy Q x ’ 

(3.18) 

and Su, S 31 are Reynolds stresses defined by 



„ dU 2 , dU x V x , dUxWx 

11 ~ dx 1 dY 1 dz ' 

(3.19) 


c dUxWx , , dwi 

31 5a; 5E dz 

(3.20) 

S 11 and S 31 can be rewritten as 


c _ c( 0 . 0 ) 

-■’ll — ‘-’ll 

+ 5j°' 2) cos 2/3 z + S[ 2 { 0) E 2 + S{l 2 ) E 2 cos 2(3 z + c.c. , 

(3.21) 

0 _ c( 0 , 2 ) 

^31 — 31 

sin 2 (3z + S< 2 ' 2) £ 2 sin 2(3z + c.c. . 

(3.22) 

After some calculation, we find that 


5 ( 0 , 0 ) 

= UaU 2 sm 2 9 , 

(3.23) 

5 ( 0 , 2 ) 

= -iaC7 y 2 sin 2 ^i*iy 0 {1) , 

(3.24) 

si?- 0) 

= -zaf ? 2 sin 2 0 [AtV o (1) + 2 sin 2 > 

z 

(3.25) 


= -iaU 2 sin 2 9AW^ ) , 

Zi 

(3.26) 

C( 0 , 2 ) 

°31 

= i/3C7 2 COS 2 0[i* W 0 (1) + 2 Sin 2 W 0) ] , 

£t 

(3.27) 

q( 2 » 2 ) 

°31 

= ^/3£/ 2 cos 2 0ifF o (1) . 

(3.28) 
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From inspection of (3.21), (3.22) and the right-hand side of (3.17), we conclude that 
V 2 has a solution of the form 

V 2 = F 2 (1) £ cos (3z + t> 2 (0,2) COS 2 f3z + V$ 2 ' 0) E 2 + c.c. . (3.29) 

The fundamental component Vj is driven only by the linear forcing term, i.e. 
L \\ 1 = iaUyyA. This is exactly the same as in the two-dimensional case (see e.g. 
Wu &: Cowley, 1992). By analogy, we obtain the jump conditions 

a j ~ a J = ™>ASgn(f7 v ) , (3.30) 

dj - dj = —irirjbjSgn(U y ) . (3.31) 

Substituting (3.23)-(3.28) into (3.17), we find that the functions V^ 0,2 ^ and V^ 2,0 ^ 
satisfy 


{ Jr " = -i5 3 sin 2 6?[i*W^ 2) + 4sin 2 fW„* (0) IF 0 (1) ] , (3.32) 

L^V&y = iS 3 sin 2 d\AW^ + 4 sin 2 0PF o (O) Wo (1) ] , (3.33) 

respectively, where we have put 

5 = OtUy. 

Inserting (3.14) into (3.32) and (3.33), and integrating these equations, we find the 
solutions 

A ( , _ /»-foo f + OO 

^ 2 ,yy — ~iS 3 sin 2 0 j / ^ 0,2) (^^M*(<i-*/M(<i-»y-Oe“ <n ^»7 , 

J 0 v 0 

(3.34) 

A , , _ /*H~oo r + oo 

V^yy =t'5 3 sin 2 ^ J o I?'%,v)Mh-v)Mti-l-t)e- m+2 Mdr) , 

where we have defined 

Ii°- 2) (Z,V) = + 4sin 2 0 f e~ 2,l ^~ 

Jo 

= / (, »((.7)r+4 s in , e/’’« + 2( 

Jo 

I { °\(, 7/) = e -*l« , + 3 ^) f 

= e _,1 ^ 1+3 ^ 2,7+6 ^ T,^+4,,^ ) 

In simplifying V^yy into the present compact form, we have dropped a purely imag- 
inary part; from (3.29) this does not alter the physical velocity. A similar procedure 
will be followed on solving for fF( 0 ' 2 ) later on. 


(3.35) 

■ 3<l « 2 dC] , (3.36) 

;) e 2 ..C’+3 , (3.37) 


(3.38) 

(3.39) 
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It follows from (3.34) that as Y — > ±oo 

(t> 2 (0l2) + c.c) ~ ±S5sm 4 dirj+ 00 £e- 2 -'< i \A{t 1 -Ti)\‘ 2 d(d V 
+ {order-one ‘no-jump’ terms} + o(l) . 

We conclude that a longitudinal vortex component must be included in the O(eJ) 
outer expansion of v so that a match with V^ 0 ' 2 ^ can be achieved. 

The function U 2 satisfies 


L 0 U 2 ,y = -UyyVx + F X {Y) 


dWi 

dz 


d 2 Ui 




sw 2 


dz 


where 

and 


F 1 {Y) = UyyY + U Y rTl , 


F,(Y) = + U yr r t Y + ^U tt t; . 


(3.40) 

(3.41) 

(3.42) 


The solution has the form 

U 2 = U^E cos /3z+Ui 0fi) +Ui 0,2) cos2(3z + U? fi) E 2 + U? >2) E 2 cos2/3z+c.c . . (3.43) 

In order to derive the amplitude equation, we need only the mean-flow distortions 
[> 2 (0 ' 0) and C/ 2 0,2 \ These satisfy the following equations respectively 


= -\a-'S 3 sin' 6 A-wP , 

r ^ \ ® 1 7y(°. 2 ) _ g(°.a) _ it y(9i 2 ) 

[ dh x dY* lUt - y "ar 11 y 2 ' Y ’ 

where 5{j’ 2 ^ is defined by (3.24). The solutions are 


(3.44) 

(3.45) 


1 /* OO /* OO 

Oft 0 ' = / ei m (i,n)A\h-v)A(h-r,-i)e- m, -did V , 

* 2 */o Jo 


(3.46) 


1 y» -f- QO J* “f- OO 

(/<«> = -V'S’sin’f)/ / /<»■»)«,, )^(J,-,)yl(i,-ii-0«- i0 ‘<i{*( . 

* 2 Jo Jo 


(3.47) 


where we have put 

ll°' 2 \t,T,) = I l0) (t,v)[? + 2^ + 4sin 2 ^ f 2 ( v - C)e- 2 ‘^- 3 ^ 2 dC] . (3.48) 

* 0 

Integrating U^y^ once with respect to Y , we find that as Y — > ±oo 

(f7 2 (0 ' 2) + c.c) ~ T 8a- 1 5 2 sin^7r/ +0 °r(7 7 - C)e" 2 '‘ c * | A(i x - i?)| 2 d<<fy 

J o 2 0 

+{order-one ‘no-jump’ terms} + o(l) . (3.49) 
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This is in fact a streamwise slip velocity generated by nonlinear interactions in- 
side the critical layers. The outer expansion for u must match to this; hence the 
leading-order outer expansion for u, (2.24), must contain a component representing 
a spanwise-dependent mean flow. 

The spanwise velocity W 2 satisfies the following equation 

L Q W 2 = - F 2 (Y)^ - S 31 , (3.50) 

where F 2 is defined by (3.42). The solution W 2 has the form 

W 2 = W^E sin (3z + W$°' 2) sin 2 (3z + W^ 2 ' 2 ) E 2 sin 2 (3z + c.c. . (3.51) 

Since at next order we are only interested in the interaction that generates the 
fundamental, W ^ does not need to be worked out. Moreover, we find that as far 
as deriving the amplitude equation is concerned, it is sufficient to solve just for the 
mean component W^ 0,2 ^. It follows from the continuity equation that 

W(°'V = -(2 (3)- l V$y ] , (3.52) 

and hence from (3.52) and (3.34) that as Y — > ±oo 

* -J- OO A T) 

(W {0 ’ 2) + c.c) ~4sm A e j3- l Y~ 2 I [ e~ 2s ' <3 |A(«i - v)\ 2 d(dTj + o{Y~ 2 ) . (3.53) 

io io 

An important point to note is that all the spanwise velocity components generated 
by the nonlinear interactions are bounded as Y — > ±oo. This is in contrast to the 
case of a purely viscous critical layer where unbounded growth of these components 
can occur (see e.g. Hall & Smith, 1991). This difference arises because the inclusion 
of unsteadiness in our critical layers means that as Y — > ±oo the balance is between 
the unsteady inertial term and the nonlinear forcing terms. In a steady viscous 
critical layer the balance is between the viscous diffusion term and the nonlinear 
forcing terms. For instance, in the case of W^ 0>2 ^ the balance is between Wyy^ (in 
our notation) and a nonlinear forcing which decays like Y ~ 2 as Y — > ±oo. Therefore 
grows like ‘Y + log Y\ In a wave- vortex interaction, this unbounded growth 
is one of the reasons why small- amplitude three-dimensional disturbances are able 
to generate order-one mean-flow distortions (Hall & Smith, 1991). 

The unbounded growth of a mean-flow distortion away from a viscous critical 
layer has been noted in other stability problems, e.g. for two-dimensional distur- 
bances in stratified shear flows (Churilov & Shukhman, 1987b). However, often the 
evolution of the disturbance is sufficiently rapid that a linear diffusion layer is estab- 
lished between the critical layer and the outer region to eliminates this growth (see 
also Brown & Stewartson (1978), Haynes & Cowley (1986) and §4.4). In the nonlin- 
ear wave— vortex interaction (Rayleigh) model of Hall & Smith (1991), the diffusion 
layer merges with the outer region because the amplitude evolution is sufficiently 
slow - in our notation Hall k Smith (1991) effectively have two time scales, namely 
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the time t to describe the rapid oscillation of the disturbance, and the time r to 
describe both the evolution of the basic flow and the growth of the disturbance. 
The relationship of the present work to the wave-vortex interaction approach will 
be examined in §4.4. 

We now proceed to derive the amplitude equation. For this purpose, it is suffi- 
cient to seek the solution for V 3 only. This term satisfies 


where 


LoV 3iY y = L 1 V 2 + L 2 V 1 + —[—S 12 + q- z S 32 } + ... , (3.54) 


it = -llO^ + U^nY^U^Y + ^UrrrT^-g^ 

Q 

“1“ U yyr T\ ] 

- { w t +{u ‘ Y+UTTi) §i ]i ]Li + j& ] ’ (3 - 55) 


the Reynolds stresses S 12 and 532 are defined by 


r\ o r\ 

Si 2 = ^(2 UiU 2 )+Qy{UiV 2 + U 2 V 1 ) + —(U 1 W 2 + U 2 W 1 ), (3.56) 

o o A 

S 32 = ^(UiW 2 + U 2 W 1 ) + —(W 1 V 2 + W 2 V 1 )+—(2W 1 W 2 ) , (3.57) 


and for brevity only those terms that contribute to the jumps across the critical 
layers are explicitly included. 

At this order it is only necessary to find the fundamental component, i.e. 


F 3 = V 3 E cos (3z + c.c. + . . . . 


Similarly we write 


* c , & 

"o ^ 12 J 

OX OZ 


ME cos j3z + c.c. + 


5 


with the non-fundamental components being omitted. 

We note that the relevant part of the linear forcing term, i.e. (Z iV 2 + L 2 Vi), is 
the same as F^(Y,ti)E in Wu & Cowley (1992). Thus the solution forced by it, 
denoted here by V 3 E , has the same asymptotic behaviour 





(“jPj + 2 qjbj + | p)bj)Y + (ajrj + pjdj + ^Z^log |F| 

+{±^7riSng(i/ y )(a+r j + Pj df + sjbj) + ...}. 


(3.58) 


It proves convenient to write the nonlinear forcing term M as a sum of the four 
terms, namely 

M — Mi + M 2 4" M$ + Mq , 
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where 


M 1 

= iaAU$' 0) + ^iaAU$ 2) + l -pAW* 2 $' 2) , 

(3.59) 

m 2 

= iaU 1>Y V 2 * {0 ' 2) - 2a 2 U 1 U^°’ 2) , 

(3.60) 

m 3 

= 2 iaU* tY V 2 {2 ‘ 0) , 

(3.61) 

Mo 

= iaAU { 2 °Y 0) + \iaAUi°y 2) + IpAwffi + iai*U 2 (2,0) 

Z z 



+ \aA*Ui 2 ' 2) + \pA*W?‘ 2) + iaU liY V}°- 2) - 2 a 2 U 1 U 2 (0,2) . 
z z 

(3.62) 


The nonlinear forcing terms are split-up in this way in order to aid the calculation 
of the solutions and their asymptotic behaviour. 

We now let denote the solution driven by Mj {j = 0, 1,2,3). Then 


where 

L { o ] V 3 % = M jiY , (3.64) 

and is defined by (3.12). In addition we denote the Fourier transform of V 3 * YY 
by Tjik): 

/ + OO A 

V 3 , YY e~ ikY dY . 

-OO 

To obtain the jump we need to evaluate {V 3 y(+oo) — V 3 y{~ o°)}, or equivalently 

r;(0) = V$(+°°) - ^ 3 ,y(-oo) • (3.65) 

A calculation shows that V 3 0 ^ makes no contribution to the jump. 

The forcing M\ represents the Reynolds stress generated by the interaction be- 
tween the vertical velocity of the fundamental wave and the induced spanwise flow. 
Using (3.65), and solving (3.64) with j = 1, we find that 

V$(+<») - V$(-oo) 

/*+oo/*-f-oo _ % 

= S*j 0 sui 2 e[ I , (3.66) 

Jo Jo 

where we have put 

jo = ir|5|“‘ , 

£<»«,!)) = X (0) «,i)){2{ 3 + f ! 7 +2sin ! « I\e + , (3.67) 

J 0 

and 


J?( 0 )(£,j?) = e -*i( 2 U+3( ,, j) 


(3.68) 
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The forcing M 2 represents the interaction between the streamwise velocity of the 
fundamental wave and the induced spanwise flow. In order to overcome a technical 
difficulty in evaluating the asymptotic form of V 3 y we write 


V&Y = Q 2 (y,tl) - iocU^lhyyU;^ . 
Substituting this and (3.60) into (3.64), we find that Q 2 {Y y t x ) satisfies 

Li 1) Q 2 = N 2 (Y,t 1 ) , 

where 


(3.69) 


(3.70) 


Since 


N 2 = M 2 ,y + iocV; x L%Xih t YYU?' 2) ) ■ 


L { o ] ihy = —iocUyUi , 


(3.71) 

(3.72) 


we obtain, by differentiating with respect to Y , that 

L^Ui.yy = - 2 iatJylhy • 

From use of (3.72) and the complex conjugate of (3.45), we find that 

N 2 = iaU liY V$’ 2) -2a 2 U x f)*? ' 2 ) -2zAa[7 y 1 ^ ljYyy U^ YY S’i*i° ,2) . (3.73) 
Then after solving for Q 2 {Y,ti) from (3.70) using Fourier transforms, we obtain 

<$(+<*>) - <$(-<*>) 

y»+oo/* + oo A 

= 25 i j 0 sm 4 e / , (3.74) 

Jo Jo 

where we have put 

= |x<°>((.^) J o t m<i+KH((+2’i+2()h- 3 - M ‘dC 


+ 


and 


= e“’ l[<£3+ ’ 73+(f+ ' ?)I1 . 


(3.75) 

(3.76) 
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The forcing M 3 is the Reynolds stress generated by the interaction between the 
streaxnwise velocity of the 3-D fundamental and the £ 2-D harmonic’ generated by 
the interaction of the fundamentals. On writing 

v,% = i) + v;'Kyy^T . < 3 - 77 ) 

and substituting this and (3.60) into (3.64), it can be shown that Q z satisfies 

L { J ) Q 3 = N 3 (Y,t 1 ), (3.78) 

where 

N 3 = M 3 ,Y — Uy 1 L^\U* t YY^2X°^ ■ 

Observing that 

L i o ) V 2 i y >) = 2 ,iaU y V} 2fi) + 2 iaS[ 2 { 0) , 
and using the complex conjugate of equation (3.72), we find that 

N 3 = -2iaU; l UlYYS;{ 2 ' 0) + 2XU; 1 UIyyy , (3-79) 

where is defined by (3.25). After solving for Q 3 from (3.78), it can be shown 

that 

^3,y(+°°) ~ ^j,y( — °°) 

= 45 4 josin 4 0/ + 7 +O °^ 3) U> 7 7)^( i i-O^( i i- 7 7-O^i- 2 ^- 7 ?)^ * ( 3 * 80 ) 
Jo Jo 

where 

Kl 3 \(,v) 

= K n '(i,ri) | j‘ 7)C(1 + 2AS ! (f + ») + ()‘(( ~ ()] n 0 ((,r,,()d( 

+ 4 sin' ej‘d( n 0 !?,r,,C)^(C-u)il+2A.S , (f+r|+<) ! ({-C)]e , ' (:! ”’ +3 "” ,| *,j , 

(3.81) 

and 

U 0 (Z,V,O = e - Jl(4<5+6 ^ 2+9 ' ,<J+6 ^ <+6 ’ j2<) ■ (3.82) 

By matching V z< y with the outer expansion we find that 

c+ - cj = V 3 ,y(+oo) - V 3 ,y(-oo) . (3.83) 

Combining (3.64), (3.58), (3.66), (3.74), (3.80), together with (3.67), (3.75) and 
(3.81), we conclude that 

c f — c~ 
c ] : 

= 7rfSgn(l7 y )(at rj + pjdf + Sjbj) 

*-j-00 /* -f- oo A 

+7ra 3 |{7 y | 3 sin 2 #/ [ Kj(^.,T}\X)A(ti — ()A(ti—^ — T})A*(ti — 2^—T])d^di] , 

Jo Jo 

(3.84) 
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where 


m.vw 

= kW((, V )( 2( 3 + ev) 

-t-2sin 3 0 { jr< 0 >(f,7 ,) j’\e + 2{(, - C)]e- 2, ‘ ( ’- 3 ' ,K ’<i( 

+K W ((,v) / < [C(27,+3C)--J(i+2 V +2C)]e- 3 *‘ <ca <iC 

Jo 

+2^ 1 )(^,J?) f 1/C[1 + 65 1 (^-C)(^ + r/ + 0 2 ]n 0 (^,T/,C)^C 

Jo 

+ K<'X(,v) / o i [('7+()(i!+3C)-(«-H;)(f-W-2C)]e- 3 -<^« J «l<</cJ 

l) j (v-i)-C)[l + 6 »,({ - C)CV‘ ( M<) Aj 

+2*< 1 >(f,,) Ac-^ll + 6»|« - 0(( + V + 

Jo Jo 

+ ^ (1 H6^)/^Ce- 3,,( ^ )(2 ^ K j[ C (t;-C)[l+6 5 i(e-C)(^+C) 2 ]e‘ Jlt2v,4a( ^ )vJ1 ^| . 

(3.85) 

Here the suffix j refers to the j-th critical layer, and the dependence on A is through 
Si (see (3.15)); K^\ K W and n 0 are defined by (3.68), (3.76) and (3.82) respectively. 
Although the kernel Aj(£, 77 1 A) is algebraically complicated, nevertheless it simplifies 
to the following form when A = 0: 

KtiiV) = ( 2 £ 3 + t 2 v) - 2 sin 2 #(2£ 3 - £ r f) - 4sin 4 0(^ 2 7/ + £t] 2 ) . (3.86) 

This is just the kernel obtained for the inviscid case (GC; Wu, 1991). 

We note that both (a+ — aj) and (dj — dj ) correspond to the classical ±7r phase 
shift in the outer expansion, while (c+ — cj ) is modified by nonlinearity. It is through 
this modification that nonlinear effects control the evolution of the disturbance. 


3.1 The Amplitude Evolution Equation 

By inserting the jumps (3.30), (3.31) and (3.84) into (2.35), we obtain the amplitude 
equation 

/7 A /*+oo r+oo 

— = g oTl A+ I j '}2g j K j (Z,r]\\)A(t l -t)A(t l -(-ri)A*(t l -2Z-ri)dZdT] , 
dt\ Jo Jo 

(3.87) 

where the sum is over all critical layers; e.g. for most of the right-hand branch of 
curve A of figure 1 there are two. The kernel Kj(£,7]\\) is defined by (3.85), while 

So = fo/f , 91 = fj/f , (3.88) 
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and 


fj = — 7 ra 3 sva. 2 9b 2 j\bj\ 2 \Uy 


(3.89) 


The constants / and /o are the same as in the two-dimensional case (Wu & Cowley, 
1992), namely 





(3.90) 


fo = 




2r Uyy T 


UyyUyr 


(UyUyyy ~ U*)U T 


\u y \ 2 


Uy\Uy \ 2 




+ 2 b 2UyyUr 

3 U 3 


-J‘ 


2 J 


(3.91) 


while J\ and J 2 are defined by (2.36) and (2,37) respectively. The coefficients gj 
are plotted against the wavenumber a = (a 2 + (3 2 )* in figures 3 a and 36 for the 
right-hand branch of curve A. The plot of the scaled coefficient (50 sec 0) against a 
is the same as shown in figure 3 of Wu & Cowley (1992) (provided that a there is 
replaced by a). For our purposes here it is sufficient to recall that the real part of 
g Q is always negative. 

In the inviscid limit, A — 0, the amplitude equation becomes 
dA r+°° r+°° 

^= 90 ^ + 9 / Kti^Aih-OAih-Z-tfA'it^-vWdr, , (3.92) 

ati «/o Jo 

where the kernel is defined by (3.86), 

* = £»> (3 93) 


and the sum is over all critical layers. Note that although in the present problem 
Uyy(y c ) ^ 0 , the nonlinear kernel K(^rj) is exactly the same as that of GC. This 
is because nonlinear interactions inside the critical layers are only associated with 
the pole singularity in u and w , compared with which, the logarithmic branch-point 
singularity associated with U yy (y c ) 7 ^ 0 is much weaker. In appendix B we show that 
the present results are directly applicable to the spatial evolution of a disturbance 
in a free shear layer when viscous effects are important. 

In the spirit of Stewartson &; Stuart (1972), Churilov h Shukhman (1988), GC 
and others, we require solutions to (3.87) to match with exponentially growing linear 
modes as t\ — 00 , i.e. 


A —> A 0 e 9 ° ritl as h -00 , 


(3.94) 


Following GC, the parameters Ao and T\ can be scaled out by introducing the 
rescaled variables 


A 

= Ae t( ' To+90iTltl) \g\ l * /{g 0r Ti) 3 , 

(3.95) 

t 

— gQrl~lti — t 0 , 

(3.96) 

A 

= A/(sorTl ) 3 , 

(3.97) 


22 



where g 0r and goi are the real and imaginary parts of g 0 respectively, and To and t 0 
are chosen so that: 

— * to -f ^Tq = log [^ols 1 1 2 /(^or T i ) 3 ] • (3.98) 

The evolution equation and the asymptotic behaviour then become 

at Jo Jo j 

(3.99) 

A — > e f as t —* — oo , (3.100) 

where Uj — <7j/|g|, and we have written A in the kernel as A. 

If there is only one critical layer, then \vi\ = \u\ = 1. Thus in this case the 
amplitude equation depends only on A, d and arg v. In this sense, the rescaling 
(3.95)-(3.97) achieves the same purpose as Shukhman’s (1991) introduction of a 
‘logarithmic time’. 


4 Study of the Amplitude Equation 


4.1 Finite-Time Singularity^ Structure 

A singularity structure for the amplitude equation (3.92), i.e. the inviscid limit of 
(3.87), was proposed by GC in their study of the spatial evolution of disturbances in 
a free shear layer. They showed numerically that solutions developed a singularity 
at a finite distance downstream, or in terms of our temporal evolution problem, the 
solutions blew up within a finite time. The singularity proposed was 


A ~ 


do 

w. z w* 


as 


t -> t. 


(4.1) 


where cr and a 0 are real and complex constants respectively. Although this finite- 
time singularity was identified for the inviscid case, substitution of (4.1) into (3.99) 
shows that it is unaltered at leading order by viscous effects. We obtain 

(3 + t V) = X:^Do|ao| 2 , (4.2) 

3 


f 


where 

Do = + + . (4-3) 

Jo Jo 

and |ao| and cr can be solved from (4.2). The singularity time t 9 can be determined 
numerically in the same way as described in Wu & Cowley (1992). 
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4.2 Viscous Limit 


Possibly the most surprising result of this paper is that the amplitude equation 

(3.87) , or equivalently (3.99), does not admit an equilibrium solution. This is a 
significant difference from the two-dimensional viscous case (e.g. Goldstein & Leib 
1989; Wu & Cowley, 1992). The reason for this is that the integral of the kernel 

/*+oo y + oo 

I I Kji^TjWdtdr] 

Jo Jo 

does not exist. 

In order to shed light on this observation, we now turn to examining the am- 
plitude equation under the very viscous limit A — > +oo. This corresponds to the 
situation where the viscosity is relatively large, or the disturbance is relatively weak 
(see (1.3)). In Appendix A it is shown that in this limit, the amplitude equation 

(3.87) can be reduced to 

A -3 "— - = g 0 riA + A l gA f |A(£j — v )\ 2 d 7 ] + 3 A 3 ) , (4-4) 

at i Jo 

where i\ is defined in (A. 7). The complex constant g is 

9 = y sm ! 0[l - 2sin^]r(i)i:^-j , (4.5) 

where the sum is again over all critical layers. The parameter A can be scaled from 
the leading order equation by the transformation 

A = A , and r x = X ~*T\ , (4.6) 

to obtain 

V 

-A = g 0 ? 1 A + gA[ |A(Fx — r})\ 2 dt] . (4.7) 

dti Jo 

We note that in this limit the critical layers are of width O(Aiej), i.e. 0(72“ >), and 
are thus viscous to leading order. However, sandwiching each critical layer there are 
thick diffusion layers of width ©(Ajcj ) where unsteadiness effects are still important 
at leading order (cf. Brown & Stewartson, 1978). 

A remarkable feature of (4.7) is that the nonlinear term is non-local so that the 
amplitude equation is not of Stuart-Watson-Landau type. 2 In fact the nonlinear 
effect comes only from the interaction between the fundamental and the induced 
mean-flow distortion - the higher harmonics play no role at leading order. This is 
in contrast to the case when A = G(l), where both the mean- flow distortion and the 
higher harmonics contribute. Note also that when 9 = the coefficient g is zero 

2 The authors are grateful to Prof. S.N. Brown for discussions without which the authors would 

not have realised this point. 


24 



so that the nonlinear term vanishes. The same behaviour also occurs in the inviscid 
case (GC); however, when 0 < A < oo, Kj(£ y rj\X) is non-zero even when 0 = 

Although a numerical investigation of (3.99) seems necessary, we note that (4.7) 
can solved analytically. First we scale out various constants by a transformation 
similar to (3.95)— (3.96): 


B = /(gorn) , 

^1 = gOr?lil — ^10 , 

where fro and iio are chosen so that 

— ?io + ii io = log [A 0 \g r \ 7 / {gor^i )] > 

and g r is the real part of g y i.e. 

2 

9r = y sin 2 ^[l - 2sin 2 0]r(i)i?e{]^ 5j /?7 5 } . 

The evolution equation (4.7), and the asymptotic behaviour, then become 

^ = B + |Bft - ,)!’<*) , 

dti \g r \ Jo 


(4.8) 

(4.9) 


(4.10) 


(4.11) 


and 
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(4.12) 


(4.13) 


By multiplying (4.12) by B* y and then taking the complex conjugate, it is straight- 
forward to show that 


Ml 

du 


= 2|5| 2 + 


fi ' 5 ' 2 r 

lffr| JO 


\B(ii-v)\ 2 dg 


(4.14) 


On solving a differentiated version of (4.14) we obtain 

16exp(2?i) 

( 4 — ]fct exp(2? 1 )) 2 


|J9| 2 = 


(4.15) 


From (4.15) it is clear that the sign of g r plays an important role in determining 
the terminal behaviour of solutions to (4.7). In particular, the solution develops a 
finite-time singularity when g r > 0, but decays exponentially at large times when 
g T < 0. Since (4.7) is a limiting form of (3.87), this suggests that when g T > 0, 
solutions to (3.87) or (3.99) will develop a finite-time singularity no matter how 
large A is. However, when g r < 0 it seems likely that solutions will terminate in 
a finite-time singularity if A is not too large, but will decay exponentially once A 
exceeds a critical value. Our numerical results demonstrate that this is indeed the 


case. 
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4.3 Numerical Study of the Amplitude Equation 

We have integrated the rescaled amplitude equation (3.99) using a finite- difference 
method. Two independent schemes have been used as a check: a Milne’s (predictor- 
corrector) method and an Adams-Moulton (implicit) method. Both schemes have 
sixth-order accuracy. The kernel is evaluated numerically using Simpson s 

rule. 

assume that A can be approximated by the linear solution when f ^ To? 
where To is a ‘big’ positive number; a suitable choice is determined by trial and 
error. The integral over the infinite domain D = [ 0 , +oo) x [ 0 , +oo) is approximated 
by that over a large but finite domain Do = [0 , Ao] x [0, ho]. For the viscous case, 
the tail over the domain T>i = (D - D 0 ) is neglected. This is justifiable because the 

decays exponentially as £ and/or 77 tend 
to + 00 . Different values of X 0 and Y 0 were tried in the program before deciding on 
suitably large values; we find that it is sufficient to take Xq — ho — t To. For the 
inviscid case, we approximate the tail over D x analytically using the linear solution. 
However, we find that dropping this tail has little effect on our results. 

The inviscid version of (3.99) has been studied by GC. Their numerical results 
confirm the singularity structure (4.1). Here we integrate the inviscid amplitude 
equation using the coefficients calculated for the Stokes layer. Figure 4 shows results 
for the wavenumber ck = 1 .2 on the right-hand branch of curve A. Four propagation 
angles 6 were investigated: 0 = 15°, 30°, 60° and 75°. For 6 = 15° and 30 the 
amplitudes exhibit an oscillatory behaviour, indicating a periodic energy exchange 
between the disturbance and the basic flow. Similar behaviour has been observed 
elsewhere (e.g. see GC and references therein). Local singular solutions are displayed 
as dotted lines. They show that a finite-time singularity occurs. Moreover, for 
0 - 60° and 75 °, the local singular solutions fit rather well with the corresponding 
numerical results over a substantial range of time, even though they were expected 
to be valid only near the singularity time. For the inviscid case, we have also worked 
out a power-series solution of the form a n e ^ 2n+1 ^- The recursion relation for the 
a n is given in Appendix C. Although this form of solution is valid when t — * - 00 , 
it is found to have a rather sizeable range of validity when truncated at high order, 
say 40-50th order. This provides a check on the numerical results (see also Wu, 
1991). For 6 = 60° and 70°, the power-series solutions were able to reproduce the 
first oscillatory cycle. 

We note that at Q — 60°, a resonant triad interaction can occur if the initial 
disturbance includes an O(e^) two-dimensional eigenmode with a wavenumber 2ot 
and a phase velocity c (Goldstein k Lee, 1992; Wu, 19926). However, in the present 
problem this resonance does not occur. This is because such a two-dimensional mode 
is not present in our initial disturbance, and when excited by nonlinear quadratic 
interactions, it only has a magnitude of <3(e 2 ) - weaker than the required strength 
0 (e*) for a resonant interaction. 

We now move onto the viscous case. As illustrated in figure 5 for 6 = 60° and 
d = 1.2 on the right-hand branch of curve A, we find that increasing the viscosity 
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generally delays the occurrence of the singularity. Note that g r is positive for these 
parameter values, so the singularity cannot be eliminated no matter how large A 
is. Figure 6 shows results for 6 = 30°, and the same wavenumber a = 1.2; the 
parameter g r is now negative. As can be seen, viscosity delays the time to singularity 
formation if A is not too large. However, once A exceeds a critical value 3 (between 
30.0 and 35.0 in this case), the amplitude decays exponentially, i.e. log \A\ ~ -c 0 t, 
or \A\ e -c ° r , where c 0 > 0 is a constant. This is significantly different from the 
two-dimensional case, where viscosity causes the disturbance to saturate at a finite 
amplitude. Although the waves decay exponentially, in appendix D we show that 
rather surprisingly the ‘slip’ velocity, (D.l), generated by accumulated nonlinear 
effects grows linearly with time. This is illustrated in figure 7, where we plot a 
suitably scaled streamwise velocity jump. It follows that in the outer region the 
spanwise-dependent mean flow driven by the slip velocity grows linearly with time, 
while the longitudinal vortex equilibriates at a finite amplitude. We conclude that 
the distortion has stabilised the basic flow, with the result that the Rayleigh modes 
start to decay. We also note that the growing velocity jump can be interpreted 
as a growing ‘vortex sheet’. This is intriguing since the development of intense 
shear layers is one of the characteristic precursors to the formation of small scale 
turbulence. This in turn raises the question of the stability of the thin shear layer 
to secondary disturbances. 

Figure 8 displays results for 9 = 60° and a = 0.8 on the right-hand branch of 
curve A. For this case g T is again negative, and viscosity plays a similar role as in 
figure 6. It is worth noting that at moderate values of A, viscosity can induce rather 
violent oscillations. However, the oscillations gradually disappear as A increases, 
and ultimately the disturbance decays when A is sufficiently large. The calculations 
presented here show that both A and the sign of g T determine the terminal form of the 
solution to (3.99). This conclusion is supported by other numerical calculations using 
artificial coefficients which we do not report here. However, it is worth observing 

that because g T depends on 0, for any given wavenumber a and sign of ]T Ujj3j 1 , it is 
always possible to find some 9 such that g r > 0, i.e. such that a finite-time singularity 
can occur. In this sense, blow-up is more common than in the two-dimensional case 
(cf. Wu & Cowley, 1992). 

Of course, once the singularity occurs, our theory ceases to be valid. Neverthe- 
less, the finite-time singularity does indicate that an explosive growth is induced by 
nonlinear effects, and we suggest that this nonlinear blow-up may be related to the 
bursting phenomena observed in experiments (e.g. Merkli & Thomann, 1975; Hino 
et al., 1976). Moreover, as GC argue, the present theory does not break down until 
the amplitude of the disturbance becomes order one, at which point the flow will be 
governed by the Euler equations. 

Finally, we examine the special case 9 = 45°, again taking a = 1.2 on the right- 
hand branch of curve A as an example. This propagation angle is special because 

3 No attempt has been made to determine the critical value precisely because it is very CPU 
intensive to integrate the amplitude equation for A close to its critical value. 
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in both the inviscid (A = 0) and viscous (A = +oo) limits, the nonlinear term in the 
amplitude equation vanishes. This does not occur, however, when 0 < A < +oo. In 
figure 9, we depict evolution curves for four values of A. The solutions clearly appear 
to develop a finite-time singularity. However, as A is increased the solutions change 
significantly. As shown in figure 10, for sufficiently large A, the solutions can evolve 
into a periodic oscillation without tending to a definite limit. In order to demonstrate 
that these solutions are not a numerical artifact, results are displayed at two different 
resolutions. We note that for A = 50.0, the amplitude exhibits a rather ‘chaotic’ 
transient state before relaxing into a periodic oscillation. This ‘chaotic’ transient 
becomes less noticeable, and ultimately disappears, as A is increased further. 

4.4 The Relationship with Wave- Vortex Interactions 

A feature of critical-layer analyses is the fact that surprisingly large mean flows 
and vortices are generated, e.g. in our analysis a spanwise-dependent mean flow is 
generated which is as large as the fundamental disturbance. A feature of the wave- 
vortex interaction of Hall &: Smith (1991) is that small amplitude disturbances can 
modify the mean flow by an order one amount. It seems natural to ask whether there 
is a link between the two theories, especially since Rayleigh wave-vortex interactions 
also involve critical layers. In particular could a disturbance consisting of a pair of 
oblique waves evolve through an ‘unsteady’ or ‘non-equilibrium’ critical-layer stage 
en route to a wave-vortex interaction? This is of course related to the questions: 

1. Can a wave-vortex interaction involving Rayleigh waves be established in the 
first place? 

2. Are the weakly nonlinear neutral waves in such a wave-vortex interaction 
stable? 

Before partially addressing at least the first of these questions it is instructive 
to determine the range of validity of (4.7). From (1.3) and (4.6) it follows that the 
time scale over which the growth-rate evolves, as specified by 7i = 0(1), is 

t — t 0 ~ 0(e*72 5 ) } (4-16) 

while from (1.1) and (A. 7) the time scale over which the disturbance amplitude 
evolves, as specified by t\ — 0(1), is 

T -T Q ~0(e-iR- 1). (4.17) 

Our theory assumes that these two time scales are distinct. This is no longer the 
case if 

e ~ R < , i.e. A ~ R * , (4.18) 

since both time scales are then 0(R~*). We conclude that when 

R~ < e <C R~ l , or equivalently 1 < A « Ri , (4.19) 
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the evolution can be described by (4.7). The validity of the full integro-differential 
form of the amplitude equation (3.87) is of course restricted only by A <C R*. 

From (4.6) the disturbance amplitude corresponding to the scaling (4.18) is 
0(A*e), i.e. 0(R~*). In terms of the global Reynolds number Re = i? 2 , the distur- 
bance amplitude is 0(Re~ tt) 5 which we note is the (corrected) wave amplitude scale 
identified by Hall & Smith (1991). While a direct connection with that work is not 
possible, we note that as part of a study of weak wave-vortex interactions, Brown, 
Brown Sz Smith (1992) have examined the scaling (4.18). Although they studied an 
equivalent spatial stability problem, in terms of our notation they effectively con- 
sidered a problem with the three times scales, t = Rr , Rit and r. Their amplitude 
equation is similar to (4.7), but because the r x and t x time scales have merged, the 
coefficient of A in (4.7) becomes £ 0 i x . Also, as a result of a viscous sub-Stokes layer 
adjacent to the wall which makes an 0(R~ l i) correction to the growth rate (e.g. 
Cowley, 1987), there is an additional linear term proportional to A. 

An important property of the amplitude equation derived by Brown et ah (1992) 
is that in one limit it can be reduced to a small amplitude form of the wave-vortex 
interaction equations, while in another limit it reduces to (4.7). A mathematical link 
between the wave-vortex interaction equations, and the non-equilibrium critical- 
layer equations is thus made. 

Let us now return to the discussion started in §1 concerning a normal mode of 
amplitude e 0 introduced at a time on one of the left-hand branches of the neutral 
curves. We note: 

1. If e 0 = 0(i2~*), then the evolution of the disturbance is described by the 
nonlinear equation derived by Brown, Brown & Smith (1992); 

2. If eo = o(R~e ), but log Cq 1 Ry the nonlinear evolution of the disturbance is 
described by (3.87), or its limiting forms (4.7) and (3.92). 

In the second case we have seen either that the disturbance hits a finite-time 
singularity, or that the mean flow is stabilised with the result that the disturbance 
decays exponentially. In both cases the final behaviour is described by much more 
rapid times scales than would be necessary for the flow to evolve to a slowly-varying 
wave-vortex interaction. Thus at least for a disturbance consisting of a pair of equal 
amplitude oblique waves, it seems that for a broad range of initial amplitudes the 
flow does not develop into a wave-vortex interaction. However, note that when 
exponential decay sets in near a left-hand neutral curve, we anticipate that the 
disturbance will start to grow again over the slow timescale r x as a result of an 
increase in the linear coefficient of (4.7). The subsequent nonlinear evolution of 
such a disturbance will depend, inter alia , on the form of the critical layer when the 
flow becomes nonlinear again. 
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5 Discussion and Conclusions 


We have derived an integro-differential amplitude equation, i.e. (3.99), that describes 
the evolution of a pair of oblique waves in Stokes layers and other inviscidly unstable 
shear layers. We have assumed that nonlinear effects become important while the 
local growth rate of the disturbances is small (although not too small), and have 
extended the previous analysis of GC by including viscosity. We have also relaxed 
the condition that the critical layer(s) occur at inflection points, although our results 
remain valid if they do (e.g. as is the case for free shear layers). However, we re- 
quire that the spanwise wavenumber 0 is not too small; the nonlinearity associated 
with the logarithmic branch point singularity is then much weaker than that asso- 
ciated with the pole. A scaling argument, similar to that in Wu & Cowley (1992), 
shows that this is no longer the case, because the nonlinearity from the simple-pole 
singularity is as strong as that from the branch-point singularity, when 

0~O(et) . 

This problem has been considered by (Wu, 1992a), who allows the disturbance 

amplitude to be modulated both in the spanwise direction on a scale Z = e*z, and 

2 

in time on the scale t\ ~ e*t. 

Numerical solutions of the amplitude equation (3.99) either blow-up in a finite- 
time singularity, or decay exponentially at large times. Similar behaviour is found 
analytically for the reduced equation (4.7) that describes the amplitude evolution 
in the very viscous limit. Equilibrium solutions to the amplitude equations could 
not be found, which suggests that nonlinear effects do not render the disturbance 
neutral. 

A direct comparison with experiment is hard, since we have been unable to find 
experiments on Stokes layers studying the evolution of well— controlled initial distur- 
bances. Indeed, in most experiments the instabilities are allowed to develop from 
background noise. However, in an attempt to relate our theory to such experiments, 
let us suppose that instability modes are excited continuously. Then, because there 
are well-defined left-hand branches to the neutral curves A and B in figure 1, there 
are specific times when new modes with a particular wavenumber can be excited. 
Let us assume that these modes are excited as soon as they are ‘viable’. 

If the initial amplitudes of these disturbances are extremely small, then their 
evolution as the Stokes layer slowly changes can be fully described by linear theory 
- eventually the disturbances either equilibrate at a finite amplitude or they decay 
on reaching the right-hand branch of a neutral curve (Cowley, 1987). However, if 
their initial amplitude is slightly larger, then the evolution of the disturbances can 
become nonlinear near the right-hand branch of the neutral curve. Since the most 
rapidly growing linear modes are two-dimensional, then on the basis of linear theory, 
such disturbances are likely to have the largest amplitudes (but see Goldstein & Lee 
(1992) and Wu (19926) for the possibility of super-exponential growth of oblique 
modes when there is a resonant triad of waves). This two-dimensional case has 
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been considered by Wu h Cowley (1992). However, suppose there is a preferential 
mechanism for exciting three-dimensional disturbances, e.g. small grooves in the 
plate. Then our results show that for moderate size disturbances (a) that for a 
range of ‘obliqueness’ angles the disturbances decay as a result of viscous effects, 
but (6) for other angles a finite-time singularity arises accompanied by an explosive 
increase in wave-amplitude. For slightly larger disturbances viscous effects have no 
time to act, and for all obliqueness angles there is an explosive increase in wave- 
amplitude. 

At even larger levels of background noise a disturbance of given wavelength will 
become nonlinear at times well before the right-hand branch of the neutral curve, 
and our theory is not applicable. A further increase in the background disturbance 
level means that nonlinear effects need to be included at times near the left-hand 
branch of the neutral curve. Again, two-dimensional disturbances have the most 
rapid linear growth. However, over much of the left-hand branch of curve A, the 
two-dimensional critical layer is regular, and so nonlinear effects lead to algebraic 
rather than exponential growth (e.g. Gajjar & Smith (1985), Churilov &; Shukhman 
(1987a), Goldstein & Hultgren (1989)). Therefore, three-dimensional instabilities 
are potentially important, although a full analysis would involve testing the stability 
of the flow with a two-dimensional, quasi-equilibrium, critical layer to secondary 
three-dimensional disturbances. This is an extensive calculation, which we do not 
tackle here. Moreover, as Killworth & Mcintyre (1985) and Haynes (1985) have 
observed, the quasi-equilibrium critical layer itself may be unstable to very high 
frequency fwo-dimensional disturbances with wavelengths comparable to the thick- 
ness of the critical layer. 

Instead we note that when three-dimensional disturbances are preferentially ex- 
cited so that their evolution becomes inviscidly nonlinear near a left-hand branch, 
explosive growth of the amplitude will again occur for all pairs of oblique modes. 
However, if the initial disturbance is sufficiently large , then viscous effects mean 
that explosive growth can only occur for a range of obliqueness angles. We believe 
that this is experimentally significant because it implies that for a relatively large 
range of initial disturbances, i.e. e 0 = o(R~6 ) but log Co 1 < R * , there is a spread of 
obliqueness angles for which viscous effects cannot prevent explosive growth. This 
is so whether or not viscous effects can force two-dimensional disturbances to evolve 
into quasi-equilibrium states. For even larger initial disturbances, i.e. e 0 = 0(R~s ), 
the amplitude equation is modified to that of Brown et al. (1992). 

Our results suggest that experimental observations are likely to depend on the 
background level of disturbances. In particular, the lower the background level of 
disturbance, the later in the oscillation cycle that nonlinear disturbances should be 
observed. For example, Monkewitz Sz Bunster (1987) note that ‘the first visible finite 
amplitude disturbances appear shortly before and around flow reversal at the edge of 
the boundary layer’; our theory when applied close to the right-hand branch of curve 
A should be relevant in this case. On the other hand, Akhavan et al. (1991a, b) note 
that at relatively large Reynolds numbers ‘turbulence appeared explosively towards 
the end of the acceleration phase’. This is a little earlier that could directly be 
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explained using our theory applied to the left-hand branch of curve A, but we note 
that Akhavan et al. (1991) conducted their experiments on finite-width Stokes 
layers in a pipe, and that their ensemble-averaged velocity profiles differ from the 
laminar profiles that we have assumed. 

As in GC, an important feature of our results is that nonlinear interactions inside 
the critical layers generate in the main part of the flow both a spanwise-dependent 
mean flow of the same size as the fundamental wave, and a longitudinal vortex. We 
note that a relatively strong spanwise-dependent mean flow was observed by Hino 
et al. (1983) in a finite Stokes layer between two plates. Hino et al. (1983) also 
suggested that the vortex structures and the bursting processes that they observed 
in Stokes layers are similar to those observed in turbulent boundary layers. This is 
possibly not surprising given (a) that our theory is applicable to any high-Reynolds- 
number shear flow that supports Rayleigh waves, and (6) that the large scale coherent 
structures in the outer region of a turbulent boundary layer may generate such 
shear flows in the wall region. Thus the results obtained here may be applicable 
to the understanding of wall-layer phenomena in turbulent boundary layers such as 
‘streaky’ vortices and high-frequency bursts. 

In addition we recall that the formation of A-vortices in boundary-layer tran- 
sition has been linked with the secondary instability of Tollmien-Schlichting waves 
and longitudinal vortices to high-frequency (Rayleigh) modes, e.g. Betchov (1960), 
Greenspan k Benney (1963). Since the Tollmien-Schlichting waves and longitudi- 
nal vortices are quasi-two-dimensional in the sense that their wavelengths are much 
larger than the thickness of the boundary layer, a slight modification to our analy- 
sis should yield an alternative nonlinear approach to describing the ‘spike’ stage of 
transition. Moreover, the generalisation of the analysis to three-dimensional basic 
states should also make it possible to extend Hall k Horseman’s (1991) linear sec- 
ondary instability analysis of Gortler vortices into the nonlinear regime. We also 
note that it is straightforward to extend our analysis to compressible flow. In par- 
ticular, essentially the same amplitude evolution equation, i.e. (3.99), is obtained. 
Our results are especially applicable to such flows since the most rapidly growing 
linear disturbances tend to be three-dimensional. 

Another consequence of our work is its relationship to the (Rayleigh) wave-vortex 
interaction of Hall k Smith (1991). An implicit assumption in the Hall-Smith theory 
is that the weakly nonlinear Rayleigh wave, whose slow evolution forces the order 
one change in the mean-flow, is stable. This assumption seems reasonable given, 
for instance, that linearly unstable Gortler vortices can evolve through a weakly 
nonlinear stage to a strongly nonlinear stage in which the mean flow is distorted by 
an order one amount (Hall, 1991). 

Initially we anticipated that our linear disturbances would be able to evolve 
through a weakly nonlinear stage near a ‘left-hand’ neutral curve, before developing 
into a wave-vortex interaction at later times (e.g. as a result of saturating into an 
equilibrium state). However, we found that the disturbances either evolve to an 
‘Euler’ stage through a finite time singularity (GC), or they decay exponentially. 
Similar results have been obtained independently by Brown et al. (1992) for distur- 
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bances which become nonlinear in the | Ar| = 0(R ~ asymptotic regime near a 
left-hand neutral curve. 

We conclude that for a wide range of initial amplitudes, a pair of oblique modes 
do not evolve to a (Rayleigh) wave-vortex interaction, even though such a distur- 
bance seems a natural initial condition for such flows. It remains to be checked 
that this is still the case for a pair of oblique waves of unequal amplitude. In addi- 
tion, the effect of both streamwise and spanwise modulation of the waves should be 
investigated, especially on the form of the finite-time singularity. 

The authors would like to thank Prof. J.T. Stuart, Dr. M.E. Goldstein and Prof. 
S.N. Brown for helpful discussions. 
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A Appendix A 

In this appendix, we show that the integro-differential amplitude equation (3.87) 
reduces to (4.4) as A — * +oo. 

In order to obtain an asymptotic estimate of the nonlinear term, we first split it 
into the following sum: 

Jo Jo 

= j\r(°> + jV(i) + 2 sin 2 0[iV (2) + N&] + 8 sin 4 0M 4) + iV {5) , (A.l) 

where 

AT (0) = 2 [ +0 ° / + °°^ (0) U,7?)^ 3 A(t 1 -OA(ti-^-7/)>i*(t 1 -2C-7 ? )^d7 ? , (A. 2) 

Jo Jo 

n m = / + “ / + “ K (l »K,,)fV(i I -{)/l(( I -£-7M , (<.-2£-^)<i{^ , (A.3) 

Jo Jo 

n {2) = 

A(t 1 -()A(t 1 -t-T))A*(t l -2t-T ] )dZdT] , (A. 4) 

= Io + °° C" 27 ? +3 CHU+2^ 2 C)]^ 3Jia3 dc| 

A(t 1 -OA(ti-^-J7)A*(t 1 -2^-7?)d^77 , (A.5) 

N W = / 0 + 7 0 +O °{^’^ / o ^Ce- 3, ‘^7^-^)[l+6s 1 (e-C)C 2 ]e" ,l(2vW) ^} 

A(t 1 -OA(ti-^-J?)A*(ti-2^-J7)^dr/ , (A- 6) 

and iVf 5 l denotes the rest of the nonlinear term, i.e. the part associated with the 
fourth, fifth, seventh and eighth terms of the right-hand side of (3.85). 

In order to estimate these integrals as A -* oo, it is necessary to make one of the 
following change of variables: 

Substitution I: 

h = A*?! , < = , v = A - ’£ , (A. 7) 

i = A"j£ , t) = A *fj ; (A. 8) 

Substitution II: (A. 7) and 

i = A - *£ , r] = \~*fi . (A. 9) 
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The appropriate substitution must be chosen in order that the resulting integrals 
are convergent. We find that substitution II converts the nonlinear terms into a 
classical cubic form, i.e. the history effects are damped out. 

Applying substitution II to and and taking the limit A — > + oo, we 

have 

7V(°) _> \-\b 0 A\A\ 2 , (A. 10) 

and 

^ \-*>C 0 A\A\ 2 . (A. 11) 

where Bo and Co are constants defined by convergent integrals. However if we 
perform substitution II in N^\ the resultant integral diverges. Instead, we integrate 
N W by parts with respect to 77, and write it in the form: 


jVW = (-3si) K {a \t,T})T)A(ti-t)-^-[A(t 1 -£-Ti)A*(ti-2Z-T])]dZdr] 

* + 00 + OO 

-(33i)- 2 / / . 

Jo Jo UL\ 

(A. 12) 

We now perform substitution I and take the limit A — > +oo; we obtain 

where j3j = | cz 2 U 2 (yi ). Note that substitution I leaves the integral convergent by 
virtue of the exponential decay of A as fj — > -foo. 

In the case of N^ 2 \ we first make the transformation (t] — 0 — ) ► 7 / to obtain 


N {2) 


_ >2 e - fcl< »-a. 1 «»(, + o-2. I c*- 

- Jo Jo Jo * 




A(t 1 -()A(t 1 -(- v -C)A*(t l -2(- v -C)dCd(d v 

_l_ / / / [2^77]e — 2jl ^ ~ 3sii (v+O-t'iC s -3«i« 2 

Jo Jo Jo 


(A. 14) 


Performing substitutions II and I into the first and second integrals respectively, 
we find that the first term is order A“a , while the second term is order A -1 . More 
precisely, 


N«> = | A(i, - + O(A-i) . (A. 15) 

In order to estimate N ^ we first integrate by parts to obtain 

r+00 /» + 00 

N {3) = - e- 53 ^ 3 - 3j '^^ 3 A(t 1 -OA(f 1 -^-7,) J 4*(t 1 -2^-7 / )^d7 7 

Jo Jo 
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_ /■ +o y +0 V 6 ^ 3 - ^'evAih-OMti-t-vWh^i-rjWdr, 

Jo Jo 

*4-00 y+oo _ 

+ 1 l K m ((,v)J o (6*i(C)[i)< , +C , -£ , C-2f’jC- {C ! ]e- 3, '« 

A^-^A^-^-^A^ti^C-^dCd^drj . (A. 16) 


Performing substitution II in the first and third terms, we can show that these two 
terms tend to A - r£>oA|A| 2 , where Do is a constant defined by a convergent integral. 
The second term is similar to JVU); the only difference is that the exponent 5si£ 3 
now replaces 2sj£ 3 of (see (3.68) for the definition of K ^). We conclude that 
as A — » +oo, 

N& ~ O(A-J) . (A. 17) 

Finally for we first perform the transformation (tj — v) — > 77, and then we 
write it in three parts: 


JV«> = f + °°d( [< dr , [+" 

Jo Jo Jo Jr] 

(v-ai+esiit-OCMti-OAiti-z+v-vWh-it+v-vWv 

- f ( d(e- 3,li<1 /’ + °°e- 2j,f3 - 3,l€2(,,+u) - , * (2v3+3 ^ J) 

J 0 Jo Jo Jo 

[C + 6 si(£ — ()( 2 (tj + ()]A(ti — — £ — ?/ — v)A ( t\ — 2( — rj — v)dv 

- f +C °dr) /dCe- 3,lC(J [ + ™ e - 2aii3 - 3 ’ li ' {r,+v) -' l{2vZ+3iv3) 

Jo Jo Jo Jo 

T]A(ti-£)A(ti-Z—ii~-v)A t (ti—‘2£—T}-v)dv . 

(A.18) 


By use of substitution II we can show that the first and the second terms tend to 
(A-t£ 0 A|A| 2 ) as A — > +00; here E 0 is again a constant defined by a convergent 
integral. Use of substitution I in the third term shows that it contributes a leading- 
order term of order A” 1 : 

i 1 1 y + oo 4 

JV«) = -A- 1 /?; 5 — — r vd)A f |A(ii - f})\ 2 dfj + 0( A-») . (A. 19) 

18 x 25 3 Jo 

Combining (A. 15) and (A. 19), and using (A.l), (A.10), (A. 11), (A. 13) and (A. 17), 
we obtain the final estimate for N : 

N _ x-^;^sm 2 d[l - 2sin 2 ^]r(hA [ +0 °\A(ii - v)\ 2 dfj + 0(A"*) . (A.20) 

J 9 3 jo 

Using (A.20), (3.87) and (3.15) we obtain (4.4). 
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B Appendix B 

In this appendix we combine our results for the Stokes layer, with those of GC for 
a shear layer, to deduce the viscous version of the amplitude evolution equation for 
a shear layer. 

First, we introduce a viscous parameter A by scaling the local Reynolds number 
R = $ 0 A/i/ in the following way: 


R- 1 = Ae . 


(B.l) 


Then the amplitude equation is 

1 dA 1 /" + 00 /* + OO 

-— = A- -7tan 2 ^y o J o K{(,ri\\)A(x-t)A(x-(-Ti)A*(x-2Z-T})dtdTi , 

(B. 2 ) 

where the kernel K is defined by (3.85) - the suffix j has been dropped because 
there is only one critical layer in a free shear layer. The parameter s\ involved in K 
now should be defined by 


*1 = itat/a'A/l-is.CW ;] 3 . (B.3) 

Readers should consult GC for the definitions of 6, Si, U c etc.. The constants R and 
7 are defined by equations (3.70) and (3.71) of GC respectively. 
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C Appendix C 

As t — > —oo, the amplitude equation (3.99) has a solution of the following form: 

^ = ga»e< 2 "«»'. (C.l) 

n=0 

Substituting the above expression into (3.99), equating the coefficients of e ( 2n+1 ) r , 
and setting a 0 = 1, we obtain the recursion relation 

1 n / 

®n+l = XT ~~TT ^ a n-l a k a l-k i 

2( n + 1) i =0 k=0 j 

n = 0,l,2,..., (C.2) 


where 


Q j (k,l,n) = [ +0 °[ + °° |A)e- 2 ( 2n " ,+2){ - 2(n+fc - ,+1)T, ^d7 ? . (C.3) 

Jo Jo 


For simplicity, we only consider the inviscid limit A 
integrate Q(k,l,n) analytically to yield 


0. In this case, we can 


Qj(k, l, n) 


1 f 6(1 — 2 sin 2 (9) l-4sin 4 fl 

16 ( (n + fc — / + l)(2n — / + 2) 4 (n + fc — 1-\- l) 2 (2n — / + 2) 3 

2sin 2 0(l — 2 sin 2 #) 1 
(n + k — Z) 3 (2n — Z + 2) 2 J 
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D Appendix D 


In this appendix, we show that although the wave decays exponentially, the nonlin- 
early induced streamwise slip velocity ( 3 . 49 ) actually grows linearly with the time 
t\. Without loss of generality, we consider 

f + oo frj 

AU= I / (77 - C)e~ 2 ^ I A(h - v)l 2 dCdrj . (D.l) 

J 0 JO 

This is essentially the shear produced by nonlinear interactions inside the critical 
layer. By means of the transform (tj — 77) — » J7, we can write AU as the sum of three 


terms 


A U — ti Ai — A2 — A3 , 

(D.2) 

where 

Ai(ti) = 

r \a(v)\ 2 d V r\-^dc , 

J -00 Jo 

(D. 3 ) 


A2(tl) = 

J \ A (v)\ 2 dr] f 71 Ce~ 2 ,l<i d( , 

J -00 Jo 

(D. 4 ) 


A 3 (£i ) = 

[ v\A{ri)\ 2 dr) f " e~ 2 * 1 ^ d(, . 

J — 00 Jo 

(D. 5 ) 

Since the integrands in (D. 3 ) are positive, we have that 


Ai < 

/“ \A M \^r 

J- 00 Jo 

e~ 2 l '< 3 d( < [ + °° \A( v )\ 2 d v r°° e- 2 ’'<*d( . 

J - 00 Jo 

(D.6) 

Similarly we find that 




VI 

(N 

< 

^ + 00 p -f- 00 

\A( v )\ 2 d v / Ce- 2a ^dC - 

-00 Jo 

(D. 7 ) 

and 




|A,| < J 

r l i n r t i- 1 

\vA(v)\ dr) j 

— 00 J 0 

V 2 “< J dC < f + °° \ V A(rj)\ 2 d V [ + °° e-^'dC . 
J- 00 Jo 

(D.8) 


The estimates (D.6)-(D.8) mean that only the term involving Aj contributes to the 
leading order term in (D.l). Moreover, Ai(tj) is a monotonically increasing function. 
Hence after a little algebra we conclude from (D.6) that as t x — > -f 00, 


/ + 00 r + oo 

\A{r})\ 2 d V J o d( > 0 , (D.9) 


and that 

where k = Ai(+oo). 


A U —> nti , 


(D.10) 
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Figure 1 . — Sketch of the linear neutral diagram for wave- 
numbers a = (a 2 + P 2 ) 1/2 (from Cowley, 1 987). t 0 is a point 
on one of the neutral curves. In the analysis we concentrate 
on times close to t = tq + T lf where t Is the magnitude 
of the disturbance and t 1 is an order one number. 


y 



Figure 2.— Evolution stages and critical-layer structures (for the case when non- 
linearity becomes Important near the right-hand branch of a neutral curve). 
The disturbance initially grows exponentially according to linear theory. As it 
approaches the neutral time tq, the growth rate become small, and linear 
critical layer(s) emerge. When the growth rate has decreased further to 
q( € 1/3^ nonlinear interactions inside the critical layers control the evolution of 
the disturbance. The earlier linear, and the subsequent nonlinear evolution 
stages, match in the overlapping domain. 
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(b) Imaginary parts. 


Figure 3. — The real and imaginary parts of the scaled coefficients in the amplitude equation (3.87), where G j = gj/[a 4 sin 2 6) 
fl = 1 ,2), and G = + G 2 . 
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Figure 4.— Ln I A I vs. the scaled time t for a = 1 .2 and A = 0 (invlscid limit): (a) 0 = 1 5°; (b) 0 = 30°; 

(c) 6 = 60°; (d) 0 « 75°. Solid lines: numerical solutions; dotted lines: local asymptotic solutions (4.1) 


(a) (b) (c) (d) 



Figure 5 Ln I A l vs. the scaled time t for a = 1 .2 and 0 = 60°: (a) X = 0; (b) X = 1 0; (c) X = 30; (d) X = 50. 

Solid lines: numerical solutions; dotted lines: local asymptotic solutions (4.1). 
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Figure 6.— Ln I A f vs. the scaled time t for a = 1 .2 and 0 = 30°. X = 10, 30, 35, 50, 80, and 1 20. Solid 
lines: numerical solutions; dotted lines: local asymptotic solutions (4.1). 



Figure 7. — 'Vortex sheet’ development: AU vs. the scaled time t for a = 1 .2 and 0 = 30°. X = 35, 50,80, and 
120. Solid lines: numerical solutions; dotted lines: asymptotic slopes. 
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Figure 9. — Ln I A I vs. the seated time t for a * 1 .2 and 0 = 45°. X = 0.1 , 1 0, 30, and 45. 
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Figure 10.— Ln I A I vs. the scaled time t for a = 1 .2 and 0 = 45°. X = 50, 60, and 100: To display the graph 
clearly, the curves for X = 60 and 1 00 sire shifted upwards by three and five units respectively. The 
Integration time steps for X = 60 and 100 are 0.05 (solid lines), and 0.1 (dotted lines). The integration 
time steps for X = 50 are 0.025 (solid lines) and 0.05 (dotted lines). 
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